Geometric Spin Hall Effect of Light 
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We describe a novel phenomenon occurring when a polarized Gaussian beam of light is observed 
in a Cartesian reference frame whose axes are not parallel to the direction of propagation of the 
beam. Such phenomenon amounts to an intriguing spin-dependent shift of the position of the center 
of the beam, with manners akin to the spin Hall effect of light. We demonstrate that this effect is 
unavoidable when the light beam possesses a nonzero transverse angular momentum. 

PACS numbers: 



Introduction. — Angular momentum (AM) of light 
beams is a topic that has recently attracted the attention 
of many researcher, from both the classical and quantum 
optics communities [H. However, all the studies produced 
up to now deal exclusively with the longitudinal compo- 
nent J z of the AM. Surprisingly enough, not attention at 
all was devoted to the study of the transverse components 
of the AM of a light beam. 

With this manuscript we begin a systematic investi- 
gation of the properties of the components of the AM 
perpendicular to the propagation axis of the beam. As it 
will be shown below, these components may be responsi- 
ble for some intriguing and counterintuitive effects, as the 
"geometric spin Hall effect of light" (geometric SHEL, for 
short) presented herewith. SHEL is an interesting phe- 
nomenon occurring when a light beam impinges upon a 
planar interface separating two different media 0, HJ ■ In 
practice, when a linearly polarized beam of light is re- 
flected or transmitted by such interface, it splits into its 
two left/right circularly polarized components. This split 
occurs in a direction perpendicular to the plane of inci- 
dence: It has been recently observed for a beam trans- 
mitted across an air-glass interface 0] • The occurrence of 
a similar left/right shift, affecting polarized light beams 
propagating along curved trajectories, has also been pre- 
dicted and observed 0. 

In this manuscript we describe a third type of left /right 
shift that occurs when a linearly polarized Gaussian 
beam is observed by means of a detection system whose 
optical axis is not parallel to the propagation direction of 
the beam. We show that this shift is inherently connected 
with the existence of a transverse part of the angular mo- 
mentum of the beam. 

Linear and angular momentum of a light beam. — A 
monochromatic electromagnetic beam in vacuo possesses 
a time-averaged linear (p) and an angular (J) momentum 
density equal to 



where Re [E(r)e~ Mt ] and Re [B(r)e _iwt ] are the time- 
harmonic electric and magnetic fields of the beam, re- 
spectively jy]. The linear momentum density p is equal 
to 1/c 2 the Poynting vector. In the quantum theory of 
light, a polarized photon of energy hui propagating in the 
z direction has a z component of spin angular momentum 
(helicity) of ah,, where a = ±1 for a circularly polarized 
photon and a = for a linearly polarized photon 0]. 
Analogously, classical electrodynamics establishes that 
for a paraxial beam of light propagating in the z direc- 
tion, the ratio of spin angular momentum density flux 
along z to energy density flux along z is 



J* pu 
cP z 



(T 

! 



(2) 



where J z — z • J(z), P z — z • P(z), and P(z) and J(z) 
are the linear and angular momentum of the beam per 
unit length obtained by integrating p and j over the x-y 
plane, respectively [7, 8]. At any plane orthogonal to the 
z direction, the intensity of the beam I(r) = c 2 p z (r), can 
be regarded as the spatial probability distribution of the 
transverse coordinate vector rj_ = xx + yy. The mean 
value (rj_) = x(x) + y(y) with respect to the distribution 
J(r) is 



(r±> = 



J J r±p z (r) dxdy 



p z (r) dxdy 



(3) 



and it determines the centroid (or barycenter) of the 
beam. Note that the denominator of Eq. ([3]) is, by def- 
inition, equal to P z . From this fact and Eq. (|lb[) , it 
immediately follows that 



p(r)=e Re[E(r)xB*(r)] 
j(r) = r x p(r), 



(la) 
(lb) 



Jx = (y) Pz - ZPy, 

that, in the plane £ = 0, reduce to 

Jx/Pz - (V) , Jy/Pz = 



(4a) 
(4b) 



(5) 
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This remarkably simple result shows that the centroid 
of a beam with a nonzero transverse angular momentum 
per unit length Jj_ = x.J x + y J y and propagating in the 
z direction, is displaced with respect to the propagation 
axis z in a direction orthogonal to «7j_: 



'J-l«=o 



= 0. 



(6) 



Equation Q automatically furnishes a simple recipe to 
actually measure the transverse angular momentum of 
the beam, by measuring the position of its centroid in 
z = 0i. 

Paraxial beams. — We will calculate the linear and an- 
gular momentum per unit length P(z) and J(z), for an 
arbitrary paraxial beam of light. Let / = f(r) denote a 
solution of the scalar paraxial wave equation [10| 



d 2 J + d 2 f + 2ikdJ = 0, 



(7) 



where k = lu/c is the wavenumber. The electric and 
magnetic vector fields are expressible in terms of / as: 

E(r) = ilu [a/x + /3/y + i (a d x f + d y f) z] , (8a) 
B(r) = ik [-0 /x + afy-i (0 d x f - a d y f) z] , (8b) 

where u = a~k+0y is a complex unit vector perpendicular 
to z that determines the polarization of the beam, and 
= 10. Next, substituting Eqs. (JSaJ) 



u u = \a\" + 

and (|8b[) into Eqs. (fTJ) we obtain 

k Px = - Im (fd x f*) + ° Re (fd y f*) , (9a) 

kp y = - Im (fd v f*) - ^Re (fd x n , (9b) 

A%=fc|/| 2 , (9c) 

and 

kj x =ky\f\ 2 + zlm(fdyt) + zaRe(fd x f*) , (10a) 

kj y = - MJ1 2 -zlm(fd x f*) + zaRe(fd y f*) , (10b) 
fcj*= -xlm(fd y n+ylm(fd x n 

-a[xRe(fd x n+yRe(fd v f*)}, (10c) 

where cr = i(a0* — a* (3) denotes the helicity of the beam, 
and both p and j are given in unit of hk. It is easy to 
check that the constraint r • j = Q is automatically satis- 
fied by the expressions in Eq. (|T0|) . We note that, in both 
Eqs. ([9]) and (fT0|) . the polarization coefficients a and j3 
appears only in the combination a = i(a(3* — a* 0), thus 
permitting an unambiguous identification of the spin and 
orbital contributions to j 



n|. 



As we are interested to 
the linear and orbital momenta per unit length P and 
J respectively, we must integrate the expressions in Eqs. 
(IH) and (fl0|) over the x-y plane. To perform this opera- 
tion we have to choose a particular function /. However, 
some general features of P and J can be inferred with- 
out making such a choice. In fact, for any function /(r) 



that vanish sufficiently fast for |rj_| — > oo, integration by 
parts leads to to following relations: 

JRe(fd x f*) = JRe(fd y r)=0, (11a) 
f[xRe(fd x f*)+yRe(fd y r)] = -/|jf, (lib) 

where / g is a shorthand for the integration of g over the 
x-y plane. From Eqs. (flQ|) and (fTTj) it follows that 

P = Jf*(z-iXV ± )f, (12) 
J = Jf*[X*i + Tx(i-iXV ± )]f, (13) 

where we used the suggestive notation X = l/k = A/ (2n) 
@. Although similar expressions for some components 
of P and J have been given before 0, @, 0, EH, Eqs. 
(fT2|) and (fT"3| are novel in that they naturally embody 
the formal connection between paraxial optics and single- 
particle quantum mechanics. It is tempting to associate 
the integrand in Eq. (p~3|) with the total angular mo- 
mentum of a single-photon wave packet with the wave 
vector centered in k = kz and spatially localized around 
r: j = X (r x k+ crk/fe), where we have neglected the 
corrections of order 0(X /|rj_|) with respect to 1. This 
expression for j in fact coincides with the corresponding 
one derived in the context of geometrical optics 0, Q • 

Further insights can be gained from Eqs. (fT2|) and 
(fT"3| if we write, without loss of generality, the arbitrary 
function / in terms of Hermite-Gaussian solutions of the 
paraxial wave equation: / = X^ m =o fnm'>Pnm( r ), with 
f nm € C, and 



irw 2 {z)n\m\ 



a/2: 



w(z) 



\/2? 



w(z) 



x e 



-i(n+m+l) arctan(2/ L) 



(14) 



Here H n (u) is the Hermite polynomial of degree n, L = 
kw^/2 is the Raleigh range of a Gaussian beam with a 
minimum spot size wq, and w(z) = wqa/I + z 2 /L 2 is the 
spot size at distance z from the waist of the beam 0]. 
After a tedious but straightforward calculation we arrive 
at the following expression for P and J: 



P = -i 



n,m,p,q 



fnm B n p5 m q + y S n p B m q 



J =x 



n,m,p : q 



-\~izS n p 5<mq) fpq 3 



+ Z [aS n p5 mq - iy/nq 8 qi m+l5n,p+l 

+iy/mp 5 m ^ q+ i5 p ^ n+ i]} f pq , 



where we have defined 



Br, 



l/(kw ) (y/pSp. 



n+1 



/nS. 



n,p+l ) 7 



C np = (kw ) (y/pSp^n+i + Vn5 n ,p+i) /2. 



(15) 

(16) 

(17a) 
(17b) 



3 



From the expressions above we note that both P and J 
do not depend upon z. However, Eqs. have general 
validity and they must be satisfied for all values of z. 
This implies that (rj_) must be a linear function of z, as 
it can be easily seen by taking the derivative with respect 
to z of both sides of both Eqs. Igj), obtaining 



d(r ± )/dz = P ± /P z 



(18) 



This equation simply states that the centroid of the beam 
propagates along the axis z obeying the laws of geomet- 
rical optics. 

To enlighten our results, it is instructive to consider 
the specific function /(r) = /oo^oo + /oiV'oi + AoV'io 
that encompasses several interesting cases. For example, 
if /oo = 1; /oi — = fio we have the fundamental Gaus- 
sian beam with I = OAM, while for /oo = 0, /oi = 
i/\/2, fio — l/\/2, / reproduces a Laguerre-Gaussian 
beam LGj with OAM i = 10. From Eq. (JTSJ) it im- 
mediately follows 



Px = Oo Im (/oo/io) , P y = 0o Im (/o* /oi) 



(19) 



where 8q = 2/(fcii>o) is the angular spread of the beam, 
and P 2 = | /oo | 2 + 1 /oi 1 2 + 1 /io 1 2 ■ For sake of simplicity, we 
normalize P z = 1. If /oo € K, Eqs. (|19p show that either 
/io or /oi must have an imaginary part to guarantee that 
the axis of propagation of the beam is tilted with respect 
to the axis z [15[ . Similarly, from Eq. (|16[) we obtain 



J x =w Re(/ *o/io), 
J y =w Rc(/* /oi), 
J z =X<j + 2Xlm(f*J 01 ). 



(20a) 
(20b) 
(20c) 



Again, if / 00 S R, Eqs. (|20a|l and ()20b|) imply that in 
order to have a nonzero transverse orbital angular mo- 
mentum either /io or /oi must have a rea/ part. It is 
known that a superposition with real coefficients of the 
fundamental mode ipoo with either ipio or i/jqi describes 
approximatively a displaced Gaussian beam 15|. Thus, 



we have shown that a lateral displacement of a Gaussian 
beam changes its transverse OAM. On the other hand, 
it is also known that such a displacement does not affect 
the longitudinal OAM J z III 17 1, as it is confirmed by 
Eq. (|20c|) whose orbital part goes to zero when both f w 
and /oi are real numbers. However, for a pure Laguerre- 
Gaussian beam LGj one has /oo = and = 0. More- 
over, in this case 2 Im (/* /oi) = 1 and Eq. (|20c[) fur- 
nishes J z = A (a + 1), in agreement with previous calcu- 
lations [f| . 

Tilted beams. — Equation (|16p clearly shows that «7j_ 
does not depend from the spin ct for any paraxial beam. 
Another simple way to see this is to substitute Eq. Qllap 
into Eqs. (110a|) and (jlObj) . However, such conclusion is 
not consistent with the following argument: Consider a 
circularly polarized Gaussian beam ^oo that propagates 
along the axis z' tilted by an angle 8 with respect to the 



reference axis z. In the Cartesian frame K' attached 
to the beam, there is a unit of spin angular momen- 
tum directed along z'\ j z r = A ct, and = 0. Dif- 
ferently, in the Cartesian frame S attached to the refer- 
ence axis z, the angular momentum of the beam will 
have both a longitudinal and a transverse component 
equal to j z = Act cos 8, and \j±\ = A |o-|sin0, respec- 
tively. As the cross-section of the beam when seen from 
K is augmented by a factor 1/ cos#, we expect that the 
transverse angular momentum of the tilted beam will go 
like «7l ~ A a tan 8 in contrast with Eq. (fTB)) . Innocent 
as it is, this conclusion has a striking consequence. In 
fact, from Eq. ([5]) it follows that for our tilted beam 
|(rjj| 2=0 | cx A Icrltanfl which is either equal to zero or 
to A tan# when the beam is either linearly or circularly 
polarized. This means that the position of the barycenter 
of a Gaussian beam observed from a reference frame non 
collinear with the direction of propagation of the beam, 
changes according to the polarization of the beam! It is 
worth noting that since the two non-collinear axes z and 
z' defines uniquely a "plane of incidence", then Eq. ([6|) 
shows that the barycenter of the tilted beam is shifted 
in a direction orthogonal to such plane of incidence. In 
this respect, this effect resembles the so called spin Hall 
effect of light (SHEL) that was recently measured jij. 

Now we will support the freshman argument pre- 
sented above with rigorous calculations. First, we pa- 
rameterize the axis of propagation z' of the beam as 
z' = x sin 8 cos <f> + y sin 8 sin <j> + z cos 8 = B,(8, <f>)z, where 
R(8, (j>) — exp (0n • L) denotes the rotation matrix that 
connects K with K' . Here n = (z x z')/(|z x is a 
unit vector perpendicular to the plane of incidence and L 
is a set of three matrices {Li, L2, L3} whose component 
are [Li]jk — —tijk, where e^fc is the Levi-Civita tensor. 
Next, we choose / = Voo( r ') as solution of the paraxial 
wave equation in the the beam frame K' . We substitute 
-000 ( r ') into Eqs. ([9]) to arrive at the following expression 
for the energy density flux in K'\ 



p'(r') oc 



j x +y 



x' [x'z' - y'aL) + y' (y'z' + x'aL) 

(21) 



+z'(z' 2 + L 2 



Finally, we transform p'(r') to the frame K via the map 

p'(r') — ► p(r) = R{8, </>)p(R-\0, <P)r) . (22) 

The resulting expression for p(r) is very cumbersome and 
it will not be reported here. However, as J x and J v equals 
(y) and —(x) respectively, at z — only, we can special- 
ize to this case. Moreover, if we consider independently 
the two cases of horizontal {4> — 0), and vertical tilt 
(4> = it/2) of the beam, we arrive at the following simple 
expressions for the density of the energy flow p{x, y, z) in 
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the z = plane: 

p z (x,0,0) (xe~ 2e (l-<? o £tan0), (23) 
Pz (0, y, O)oce~ 2,)2 (1 + O »? tan 9), (24) 

were, for clarity, we defined the dimensionless variables 
£ = x/wq and n = y/w a . For a well collimated beam 
#o <C 1 and from Eqs. (|24j|23p it immediately follows 

(x) = - * (o-/2) tan 0, (25) 
(y) =*(cr/2)tan0, (26) 

in agreement with our qualitative argument. For arbi- 
trary values of z and the equations above generalize to 



the beam itself. The essence of the phenomenon is that 
if the beam possesses a transverse angular momentum, 
necessarily the barycenter of its intensity distribution is 
shifted with respect to its axis of propagation (axis iden- 
tified with the direction of the Poynting vector). Such 
a shift occurs in a direction orthogonal to the transverse 
part of the angular momentum of the beam. This phe- 
nomenon shares many features with the spin Hall effect of 
light, but it is not generated by the interaction between 
light and matter. For this reason it could be referred to 
as the "geometric spin Hall effect of light" . Experiments 
are in progress in our labs to confirm the existence of this 
curious phenomenon. 



(x) = — X (a/2) tan 9 sin <f> + z tan 9 cos cj>, (27a) 
(y) = X (a/2) tan 9 cos <f> + z tan 6 sin <j>. (27b) 

As expected for a tilted beam, d(rjj/dz = const., 
namely P±_ = const. This intuitive result can be verified 
by a straightforward calculation that produces 

P/P z =xtan6»cos</> + ytan0sin</> + z, (28) 

in agreement with Eqs. ([27)1 and (fT5)) . Finally, we cal- 
culate the angular momentum per unit length obtaining 

J x /P z =X(a/2)tan9co8(f), (29a) 
Jy/P z = X (er/2) tan 6 sin cj>, (29b) 
J z /P z = X (a/2) (2 - sin 2 9) sec 2 9. (29c) 

It is worth nothing that ( J 2 + J 2 + J 2 ) 1/2 / P z = 
(<r/2)(4 — sin 2 9) 1 / 2 sec 2 9, that is J is not conserved by 
the rotation. This is not strange because only the mag- 
nitude of the spin density j must be invariant under ro- 
tation. In fact, as the position vector r in Eq. (|lb[) is 
defined in the reference frame K and, therefore, is not 
subjected to the rotation, then the resultant J is frame- 
dependent and cannot be invariant. Finally, it is easy 
to check that Eqs. (H]) are automatically satisfied by the 
quantities expressed in Eqs. (|27H29[) . 

Conclusions. — We have presented a quite counterintu- 
itive and intriguing phenomenon that occurs with tilted 
beams of light: The position of the intensity barycenter 
of a tilted Gaussian beam, varies with the polarization of 
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